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We study the influence of deformations on magnetic ordering in quantum dots doped with mag- 
netic impurities. The reduction of symmetry and the associated deformation from circular to ellip- 
tical quantum confinement lead to the formation of piezomagnetic quantum dots. The strength of 
elliptical deformation can be controlled by the gate voltage to change the magnitude of magnetiza- 
tion, at a fixed number of carriers and in the absence of applied magnetic field. We reveal a reentrant 
magnetic ordering with the increase of elliptical deformation and suggest that the piezomagnetic 
quantum dots can be used as nanoscale magnetic switches. 
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Quantum dots (QDs) can be viewed as artificial atoms 
which allow for a versatile control of the number of car- 
riers, their spin, and the effects of quantum confine- 
ment P, 0, y|- The motivation to magnetically dope 
semiconductor QDs comes from the possibility for an en- 
hanced control of the magnetic ordering as compared to 
their bulk-like counterparts 0, S @, 0] as well as to di- 
rectly study the influence of quantum confinement and 
strong Coulomb interactions on magnetism. There is 
an encouraging experimental progress in Mn-doped II- VI 
and III-V QDs, from the controlled inclusion of a single 
Mn impurity to the onset of magnetization at temper- 
atures substantially higher than in similar bulk materi- 
als [1, [tj 13, HI, 0, EH • However, most of the theoretical 
studies have been limited to a small number of carriers 
(N) and Mn- impurities (N m ) and the effects of the sym- 
metry of quantum confinement and electron-electron (e- 
e) Coulomb interaction on the magnetic pha se diagram, 
remain to be understood 14, 15, lg, 12, HI- I n partic- 
ular, it would be useful to know what type of external 
influences can effectively produce magnetic ordering of 
carrier spins and magnetic impurities and how such or- 
dering would modify the electrical, optical, and transport 
properties of QDs. 

In the absence of magnetic doping, circular QDs with 
rotational symmetry of the lateral confining potential 
have degenerate energy levels and show a pronounced 
shell structure, similar to atoms and nuclei. The in- 
terplay of the Pauli principle and the reduction of e-e 
Coulomb exchange energy near degenerate levels in the 
single-particle spectrum, according to Hund's first rule, 
favors the parallel spin alignment of electrons in open 
shells 0, 0] ■ Breaking the circular symmetry by deform- 
ing the lateral confinement removes the single-particle 
level degeneracies and leads to transitions in the spin of 
electrons, referred to as piezomagnetism [Hll9j|. In verti- 
cal QDs, elliptical deformation can be made by micron- 
sized device mesas incorporating dots, where the lateral 
electrostatic confinement originates from side- wall deple- 
tion. Both the symmetry of the confinement and the 
effective QD size can be controlled by the action of a 



Schottky gate wrapped around the mesa in the vicinity 
of the QD TQ, 2(3, 21 1. In lateral QDs, the confinement 
potential is formed electrostatically by several gate volt- 
ages, altering the shape and the deformation of the con- 
fining potential while keeping N constant 2^, 23] ]. 

We focus here on the effects of elliptical deformations 
on magnetic ordering in (II,Mn)VI QDs. In particular, 
we are interested to explore the onset of a spin-polarized 
state with finite Mn-magnetization, in the absence of 
applied magnetic field, to which we refer as "ferromag- 
netic" (FM) state 14, 15, Since Mn is isoelectronic 
with group-II elements, the number of carriers can be 
changed by chemical doping or by electrostatic gates. 
We represent magnetic QD in zero magnetic field by the 
Hamiltonian H = H e + H ex + H m , describing contri- 
butions of interacting electrons, electron-Mn (e-Mn) ex- 
change, and direct Mn-Mn antiferromagnetic (AFM) cou- 
pling, respectively. The interacting electron contribution 
is given by 
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where we set Ti = 1, m* is the electron effective mass, 
r = (p,z), and U QD (r) = V QD {p) + V^ D (z) is the 
three-dimensional (3D) QD confining potential, where 
Vq D (z) — m*tt 2 z 2 /2 is the ID parabolic potential with 
the characteristic subband energy O, and 



Vq D (p) = V QD {x,y) = (l/2)m*w 2 (x 2 S - 
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is the 2D anisotropic parabolic potential [l| which de- 
scribes the lateral confinement with frequencies u) x = 
ujQy/S, Lo y = loq/\/S, where 5 = u x /ujy ^ 1 is the 
strength of elliptical deformation. We impose w 2 = uj x u y 
which conserves the area of QD with deformation. The 
last term in Eq. (fTJ) is the repulsive e-e Coulomb in- 
teraction screened by the dielectric constant e, and — e 
is electron charge. The e-Mn exchange contribution 
is H ex = -JsdJ2ijSi ■ Mrffa - where J sd is 

the exchange coupling between electron spin Si at r.; 



and impurity spin Mj at Rj. The Mn Hamiltonian 
is H m = J2i,r J u'Mi ■ Mr, where jff, is the di- 
rect Mn-Mn AFM coupling. The z-component of Mj 
is M z = —M, —M + 1, . . . , M, where we choose I as 
the quantization axis and M = 5/2 for Mn. We use a 
real space finite-temperature local spin density approx- 
imation (LSDA) and the mean filed approximation for 
Mn spins to avoid computational complexity [It], H3] of 
exact diagonalization, limited to only a very small num- 
ber of interacting electrons and magnetic impurities . 
An effective Hamiltonian describing electrons can be ob- 
tained by replacing the Mn spins, that are randomly dis- 
tributed, with a classical continuous field H° s = H e — 
'^2 i Jsdn m 2 f{M z {Yi)), where n m is the averaged density 
of Mn, and a = ±1 for spin up (|), and down (J,). Within 
the mean field approach (l7j . the magnetization density 
can be expressed as (Af z (r;)) = M Bm {Mb{vj) /k^T) 
where Bm(x) is the Brillouin function [251 ] . ks is the 
Boltzmann constant, T the absolute temperature, while 
6(r0 = -J c 4 F (Af,(r 4 )) + J sd [n^r t ) ~ n l (r l )}/2 is the 
effective field seen by the Mn. The first term in 6(rj) 
describes the mean field of the direct Mn-Mn AFM cou- 
pling [3] and n a (vi) is the spin-resolved electron density. 

Following the decomposition of 3D confining potential, 
we expand the QD wave functions in terms of its planar 
ipicr{p) and subband wave function (in typical disk- 
shaped QDs only the first subband is filled). We project 
i?| ff into a 2D Hamiltonian by integrating out £,(z). In 
LSDA we express the Kohn-Sham (KS) Hamiltonian as 
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where the Vh is the Hartree potential, Vxc ^ s Vosko- 
Wilk-Nusair spin dependent exchange-correlation poten- 
tial [26[ , and the exchange spin splitting is given by 
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where J ern — J s dn m M is the e-Mn exchange coupling. 
We solve self-consistently the Kohn-Sham equations, 
H K s^ia{p) = tiaipiaip), where H KS is given by Eq. © 
and ti G are the KS energies. Our numerical results 
are illustrated for parameters based on (Cd,Mn)Tc QDs 
with J sd = 0.015 eV nm 3 , m* = 0.106, e = 10.6 0, 
a* B = 5.29 nm, and Ry* = 12.8 meV are the effective 
Bohr radius and Rydberg energy. We choose a QD with 
lu q = 25.6 meV, n m = 0.1 mxr 3 , j£f = 0.02 meV, and 
having the perpendicular (z) width 1 nm. 

We first examine the effects of elliptical deformation, 
described with the anisotropic 2D lateral confinement 
from Eq. ([2]), on the single-particle levels. In the ab- 
sence of e-e Coulomb interaction and magnetic impuri- 
ties Jem — n m — 0, the corresponding single-particle 
spectrum E nx>ny (5) = uj [(n x + 1/2)VS + (n y + 1/2)/ VS], 
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FIG. 1: Effects of elliptical deformation. Top: Single-particle 
levels of the anisotropic QD as a function of deformation 
S. Degeneracies lead to the formation of closed shells with 
N = 2,6, 12, 20, .., in the isotropic (circular) case (5 = 1), and 
the formation of subshells at 8 = q/p for integer q,p. Broken 
square and the corresponding arrows schematically illustrate 
how the spin alignment for N = 6 would be modified near 
the degeneracy point at 5 — 2,E/ljo = 2.5 in the presence 
of e-e Coulomb interaction and magnetic impurities. Bottom: 
The inset illustrates the magnetic phase diagram as a func- 
tion of quasi-particle gap A 9P , and e-Mn exchange coupling 
Jem- The line A qp = J em separates ferromagnetic (FM) and 
antiferromagnetic (AFM) ordering. The magnetic phase dia- 
gram of QD with N = 6 interacting electrons at T = 0.5 K 
is calculated as a function of 5 and Jem/^o- The dashed line 
corresponds to n m =0.1 nm -3 and u>o = 25.6 meV. 



deformation 5 in Fig. [TJ top. Isotropic QD (6 = 1) have 
a pronounced shell structure with n x + n y + 1-fold degen- 
eracy in each shell, having TV = 2,6, 12,20, .. electrons. 
While, any deviation from the isotropic parabolic con- 
finement breaks the geometrical symmetry and removes 
these degeneracies, additional accidental degeneracies, 
can occur in deformed QD at 8 = p/q for integer p, q and 
lead to the formation of subshells. For example, at 6 = 2 



3 




N 



FIG. 2: The z-component of the total spin of electrons s z 
(a) and averaged magnetization per unit area {M z ) (b) as 
a function of number of electrons N at T = 0.5 K and the 
Mn-density n m = 0.1 nm -3 . Isotropic (5 — 1, filled circles) 
and elliptical QD (5 — 1.9, empty circles) show qualitatively 
different dependence of s z and {M z ) with N. 



both p- and d-, as well as d- and /-levels cross, as shown 
in Fig. [TJ top. The 5 = 2 shell structure is significantly 
modified from the isotropic case, N — 2,4,8, 12, cor- 
respond to the closed shells and N — 6, 10, 15, to the 
half- filled shells. 

For the study of magnetic ordering in (II,Mn)VI QDs, 
it is important to understand how this single-particle 
picture is modified in the presence of e-e Coulomb in- 
teraction and e-Mn exchange coupling J em . The cor- 
responding sequence of filling electrons and their spin 
alignment, through carrier-mediated magnetism, affects 
the magnetic ordering of Mn-impurities and the possibil- 
ity to externally control Mn-magnetization. 

A simple sketch of anticipated modification near 5 = 2 
for the N = 6 state, shows that the parallel spin align- 
ment of the valence electrons is favored for a wide range of 
deformation. Additional analysis is addressed in Fig. [1] 
bottom. The inset shows that the occurrence of FM or 
AFM state (having finite and vanishing magnetization, 
respectively), depends on the competition between the 
two characteristic e nerg y scales: the quasi-particle energy 
gap A qp , and J em [27j. The single-particle gap should 
be distinguished from A qp which is the energy difference 
between highest occupied and lowest unoccupied KS lev- 
els and depends on N, the total spin of electrons, the 
strength of e-e Coulomb interaction, and the underlying 
symmetries of QD confining potential. For J em = 0, as 
shown in Fig. [TJ the FM state is stable only within a 
narrow range of 5 which becomes vanishingly small for 
noninteracting electrons (at 5 — 2, Fig. [TJ top). The in- 
crease in J em extends the deformation range over which 
the FM state is possible. 

To further examine the effect of elliptical deforma- 



tion on magnetic ordering we introduce the z-component 
of the total spin of electrons, s z = (iVj — iVjJ/2, and 
the spatially-averaged Mn-magnetization per unit area 
A, (M z ) = {I /A) Jd 2 p{M z (p}), obtained from the self- 
consistent solution of Kohn-Sham equations. In Fig.[2]we 
show s z and (M z ) as a function of N for both isotropic 
QD with circular symmetry (5 = 1), and elliptical QD 
with 5 = 1.9. For interacting electrons and finite J em we 
observe that carrier-mediated magnetism leads to cou- 
pled ordering of electron and Mn spins. The oscillations 
of s z and (M z ) reveal that the shell structure is strongly 
modified with deformation. The vanishing of s z and (M z ) 
correspond to the closed shells (N = 2, 6, 12, at 5 = 1 
and N = 2,4,8, 12, at 5 = 1.9) while their local max- 
ima correspond to the half-filled shells (N = 1, 4, 9, at 
5 = 1 and N = 1,3,6,10,..., at 5 = 1.9). 

Focusing on the N = 6 state, we see that the min- 
ima in s z and (M z ) at 5 = 1 become maxima at the 
deformation 5 = 1.9. For an isotropic QD the p- levels 
are degenerate (Fig. [TJ and N = 6 corresponds to the 
AFM state with A qp > J em and vanishing (M z ). El- 
liptical deformation breaks the circular symmetry of the 
lateral confinement and removes the p-level degeneracy. 
At 5 = 1.9, J em > A qp stabilizes the FM state with 
finite (M z ) and s z = 1 implies that the electron spin 
polarization is P = (JVf - N|)/(iV T + iVj_) = (4 - 2)/6. 
The N = 6 state illustrates an intriguing possibility to 
turn the Mn-magnetization on and off without changing 
the number of carriers, needed in the bulk-like magnetic 
semiconductors. However, N = 6 is not the unique state, 
a similar behavior near 5 = 1.5 is expected also for the 
N = 12 state, forming another closed shell in the limit 
of isotropic confinement. Elliptical deformations of the 
lateral confinement in Mn-doped QDs thus lead to the 
change between the AFM and FM state which can be 
controlled by gate voltage. Such structures would behave 
as piezomagnetic QDs in which the control of deforma- 
tion determines the magnetic ordering. 

To examine if the piezomagnetic QDs could function 
as nanoscale magnetization switches, we also consider the 
evolution of their magnetic ordering with T and 5, shown 
in Fig. [3] for N = 6. While (M z ) and P monotonically 
decay with T, the inset reveals a nonmonotonic behavior 
of (M z ) and P with 5. The reentrant magnetic ordering 
with the increase of elliptical deformation corresponds 
to the AFM-FM-AFM transitions, anticipated from the 
sketch in Fig.[TJ top. Since the change between vanishing 
and finite (M z ) occurs in a narrow range of 5 (near 1.7 
or 2.1, respectively) we expect that a moderate change 
in the gate voltage would suffice to implement such a 
magnetization switch [28| . With the increase of 5 we can 
also anticipate FM-AFM-FM transitions. In the limit of 
5 = 1, a QD with N = 10 forms the FM state with finite 
(M z ) and s z = 1 but with an increased deformation and 
opening up of A gp , near 5 = 1.5, the AFM state with 
(M z ) — s z = is favored. With the further increase in 
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FIG. 3: Temperature evolution of Mn-magnetization per unit 
area (M z ) and electron spin polarization P in an elliptical QD 
with 8 = 1.9 and N = 6. Inset: Evolution of (M z ) and P with 
elliptical deformation at T — 0.5 K shows a nonmonotonic be- 
havior with sharp rise and fall (near 8 = 1.7 and 8 — 1.9, re- 
spectively) corresponding to the switching between the AFM 
and FM state which could be controlled by gate voltage. 



S, we find switching back to the FM state with s z = 1, 
near 5 = 2, also expected from a single-particle picture 
(Fig. [TJ top) and the crossing of d- and /-levels at 6 = 2. 

Another potential application of piezomagnetic QDs 
could be their voltage-controlled P to enable an efficient 
and tunable spin injection in semiconductor nanostruc- 
tures [7|. While for the N = 6 state the total electron 
spin polarization is limited to P — (4 — 2)/6 = 1/3, the 
transfer of only valence electrons could lead to the injec- 
tion of completely spin-polarized carriers, in the absence 
of applied magnetic field. Analogously, such QDs could 
also be used as voltage-controlled spin filters [29( . 

In conclusion, we have investigated the effects of ellip- 
tical deformation, electron-electron Coulomb interaction, 
and the number of carriers on the magnetic ordering in 
Mn-doped QDs. We reveal additional possibilities to con- 
trol magnetism in semiconductor nanostructures which 
would not be feasible in their bulk counterparts. The 
two previous experimental results: (1) the feasibility of 
large elliptical deformations without the Mn-do ping 11911 
and (2) controlled Mn-doping in QDs i, @, 0, [13, El [lj, 
would need to be combined in order to realize the pre- 
dicted behavior of piezomagnetic QDs. While our anal- 
ysis has focused on the electrons and II- VI materials, 



much larger hole-Mn exchange coupling |25|, |30| in either 
II- VI or III-V QDs could extend the suggested control 
and manipulation of magnetic ordering to significantly 
higher temperatures and enable utilization of piezomag- 
netic QDs for high-density magnetic memory. 
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